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Results of direct numerical simulation of isotropic turbulence of surface gravity waves in the 
framework of Hamiltonian equations are presented. For the first time simultaneous formation of 
both direct and inverse cascades was observed in the framework of primordial dynamical equations. 
At the same time, strong long waves background was developed. It was shown, that obtained 
Kolmogorov spectra are very sensitive to the presence of this condensate. Such situation has to be 
typical for experimental wave tanks, flumes, and small lakes. 
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-Introduction - In this year we have a 50th anniversary 
of the famous work by Phillips 1] which was, probably, 
the first attempt to give an explanation for power-like 
spectra of surface gravity waves observed in numerous 
experiments. In recent works [2, 3] the physical expla- 
nation given by Phillips was corrected. During less than 
ten years after that the statistical theory of surface waves 
was founded: Hasselmann derived kinetic equation for 
waves Q , Zakharov created theory of wave (or weak) tur- 
bulence p, which describes solutions of this equation. 
Stationary Kolmogorov solutions of the kinetic equation 
corresponding to flux of energy from large to small scales 
(direct cascade) and flux of wave action (waves "num- 
bers") from small to large scales (inverse cascade) were 
found [1, 01 • This opened a way to creation of the ef- 
fective tool for waves forecasting. The conjectures under 
which the theory of weak turbulence was derived includes 
Gaussian statistics for waves field and resonant interac- 
tions prevalence Q. They are subject for confirmation. 
Modern numerical methods allow to perform wave field 
modeling in the framework of kinetic equation faster than 
real processes in nature. At the same time, it is impos- 
sible to create waves forecasting model based on direct 
numerical simulation of dynamic equations. Even more, 
we do not need to know velocity and elevation at ev- 
ery point of the surface. Statistics, especially mean wave 
hight and speed, this is what really matters for estima- 
tion of operational conditions of oil platforms and cargo 
ships. And this is exactly the subject of theory of weak 
turbulence. In means that the problem of confirmation 
and correction of the waves turbulence is of great practi- 
cal importance. 

Experiments in the open sea and on the Great Lakes 
gave temporal and space spectra consistent with the the- 
ory [^, U [l^ . A comprehensive review of experiments 
and comparison with the theory of the weak turbulence 
can be found in [l^, ll^ ■ Most of these experiments 
were performed with wind pumping, broad in spectrum. 
Narrow in spectrum pumping can be realized in wave 
tanks or flumes. Results obtained on such state of the 
art devices frequently contradict with predictions of the 



theory of wave turbulence. For example, in the recent ex- 
periments [3, flit observed spectra were changing slope 
with variation of steepness and pumping force. 

May be the most promising way to check conjectures of 
the waves turbulence theory is a numerical experiment. 
In the case of direct numerical simulation we have the 
highest possible control on the parameters of experiments 
and all information about the wave field. At the same 
time all this data is given at the cost of the enormous 
computational complexity. Fast growth of computational 
power and development of computational algorithms al- 
lowed direct numerical simulation of the surface gravity 
waves, starting from the simulations of the swell evolu- 
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There is a hope, that this ap- 
proach together with confirmation of conjectures of the 
weak turbulent theory will allow us to explain phenom- 
ena observed in experimental wave tanks. 

At the same time, theory of the wave turbulence is 
still under development. To close the circle the recent 
paper by Newell and Zakharov Q gave second life to 
the Phillips spectrum, although from completely differ- 
ent point of view: Phillips spectrum considered to be a 
solution which give a balance of transfer of energy due to 
nonlinear waves interaction transfer of energy and trans- 
fer due to intermittent events, like wave breaking and 
white capping. 

This Letter was inspired by several recent papers. In 
the first one 25] numerical simulation of the isotropic tur- 
bulence with observed formation of inverse cascade was 
performed in the framework of Zakharov equations 0. 
A little bit later a group of authors [26^1 during simula- 
tion of 2D hydrodynamics observed formation of large 
scale structure due to Kraichnan's inverse cascade and 
explored its influence on the system. Approximately 
at the same time state of the art surface waves wave 
flume experiment was performed [3]. Observed spectra 
differed from the theory wave turbulence. Author per- 
formed a direct numerical simulation of isotropic turbu- 
lence of surface gravity waves in the framework of Hamil- 
tonian equations. The isotropic turbulence is a classical 
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setup for turbulence investigation. In nature such such 
wave field is usually observed in the regions with large 
amount of floating broken ice. For the first time the 
formation of both direct and inverse cascades was ob- 
served in the framework of primordial dynamical equa- 
tions. At the same time, strong long waves background 
was developed. This phenomenon of "condensation" of 
waves (following analogy with Bose-Einstein condensa- 
tion in condensed matter physics) was predicted by the 
theory of weak turbulence. It was shown, that obtained 
Kolmogorov spectra are very sensitive to the presence of 
the condensate. Such situation have to be typical for 
experimental wave tanks, flumes and small lakes. Ob- 
tained results can be considered as the first observation 
of generalized Phillips spectra, introduced in 0] and ex- 
plain some deviations from the waves turbulence theory 
in recent wave tank experiments. 

-Theoretical background - We consider a potential fiow 
of ideal incompressible fluid. System is described in terms 
of weakly nonlinear equations [1, for surface elevation 
r]{r,t) and velocity potential at the surface ^{r,t) (r = 
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-[kij]k[Tjk^P] - [rjkij]A^ + F-i[7fcVfe] + Pf- 

Here dot means time-derivative, A — Laplace operator, 
A: is a linear integral operator (^k 

inverse Fourier transform, 7^ is a dissipation rate (ac- 
cording to recent work [27] it has to be included in both 
equations), which corresponds to viscosity on small scales 
and, if needed, "artificial" damping on large scales. Pf 
is the driving term which simulates pumping on large 
scales (for example, due to wind). In the fc-space sup- 
ports of 7fe and are separated by the inertial interval, 
where the Kolmogorov-type solution can be recognized. 
These equations were derived as a results of Hamiltonian 
expansion in terms of fcry. From physical point of view 
/c-operator is close to derivative, so we expand in powers 
of slope of the surface. In most of experimental obser- 
vations average slope of the open sea surface is of the 
order of 0.1, so such expansion is very reasonable. 

In the case of statistical description of the wave field, 
Hasselmann kinetic equation Q for the distribution of 
the wave action n{k,t) = (|a7r(i)p) is used. Here 
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are complex normal variables. For gravity waves ujk — 

Vgk- 

From the theory of weak turbulence besides 
equipartions (Rayleigh- Jeans) spectrum, we know two 



stationary solutions 0] of the kinetic equation in the 
case of four- waves interaction: 



.(1) 



= CiP^/^k-^- 
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For surface gravity waves, a coefficient of homogeneity of 
nonlinear interaction matrix element /3 = 3, the power of 
dispersion law a = 1/2, and the dimension of the surface 
d = 2. As a result we get 



,(1) 



The first solution n'j^' describes direct cascade of energy 

from large pumping to small dissipative scales. The sec- 

(2) 

ond solution nj. describes inverse cascade of action (or 
"number" of waves) from small pumping to larger scales. 

-Numerical simulation - We simulated primordial dy- 
namical equations (H]) in a periodic spatial domain 27r x 
27r. Main part of the simulations was performed on a 
grid consisting of 1024 x 1024 knots. Also we performed 
long time simulation on the grid 256 x 256. The used 
numerical code was verified in [H, H^, [lH, [H [H, [ll]. 
Gravity acceleration was g = 1. Pseudo- viscous damping 
coefficient had the following form 



7fc 



0,k<kd, 

-70 (fc - kd)'^,k > kd, 
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where kj, = 256 and 70,1024 = 2.7 x 10** for the grid 
1024 X 1024 and k^ = 64 and 70,256 = 2.4 x 10^ for the 
smaller grid 256 x 256. Pumping was an isotropic driving 
force narrow in wavenumbers space with random phase: 



4P( 
0- 



(fc — fcpi)(fcp2 — fc) 
(fcp2-fcpl)^ ' 

if fc < kpi or fc > fc. 
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here fcpi = 28, fcp2 = 32 and Fq = 1.5 x 10"^ Rj:{t) 
was uniformly distributed random number in the in- 
terval (0, 27r] for each fc and t. Initial condition was 
low amplitude noise in all harmonics. Time steps were 
Atio24 = 6.7 X 10"'' and A<256 = 5.0 x 10-^. We used 
Fourier series in the following form: 
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(27r)2 



Ny/2 



-~N^/2-Ny/2 

here , Ny — are number of Fourier modes in x and y 
directions. 

As a results of simulation we observed formation of 
both direct and inverse cascades (Fig. [TJ solid line), 
although exponents of power-like spectra were differ- 
ent from weak turbulent solutions What is im- 
portant, development of inverse cascade spectrum was 
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arrested by discreteness of wavenumbers grid in agree- 
ment with 1^, 2^ 2^. After that large scale condensate 
started to form. As one can see, value of wave action |a/cp 
at the condensate region is more than order of magnitude 
larger than for closest harmonic of inverse cascade. Dy- 
namics of large scales became extremely slow after this 
point. We managed to achieve downshift of condensate 
peak for one step of wavenumbers grid during long time 
simulation on a small grid 256 x 256 (Fig. [1] line with 
long dashes). As one can see we observed elongation 
of inverse cascade interval without significant change of 
the slope. Unfortunately, inertial interval for inverse cas- 
cade is too short to exclude possible influence of pump- 
ing and condensate. We can try to estimate exponent 
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FIG. 1: Spectra < \ak\ >. With condensate on the 1024 x 
1024 grid (solid); on the 256 x 256 grid with more developed 
condensate (long dashes); without condensate on the 1024 x 
1024 grid (short dashes). 

by compensation of the spectra in a double logarithmic 
scale (Fig. The observed spectrum ~ fc~^-^ is close 
to weak turbulent solution ([3]). Slightly lower exponent 
could be explained by weakening of resonant nonlinear 
interactions on the rough wavenumbers grid, which ef- 
fectively decreases the homogeneity coefficient f3 in the 
expression ([3]). For direct cascade spectra we also used 
compensation in double logarithmic scale. Results are 
present in Fig. [3] (left). Formally, in this case we have 
quite long inertial interval 32 < A: < 256, but in reality 
damping has an influence on the spectrum approximately 
up to fc ~ 180. Still in this case we have more than half 
of a decade. Theory of weak turbulence gives us depen- 
dence ~ k~'^ ([3]), known as Kolmogorov-Zakharov spec- 
trum. Nevertheless, one can see that we observe fc~^/^, 
known as Phillips [ll, Q spectrum. So we need to un- 
derstand, what is the reason of different spectrum slope? 
What changes weak turbulent theory in this case? 
To answer these questions let us compare our situ- 



ation with previous works on decaying 17|, l2d . |21| or 



isotropic [22, 2^, 24 1 turbulence. Immediately we have an 



FIG. 2: Compensated inverse cascade spectra C{\akf)k" . 



answer: condensate and inverse cascade spectrum! The 
inverse cascade's part of the spectrum is described by 
the theory of weak turbulence, so let us concentrate on 
the strong long (fc ~ 5) waves' influence on much shorter 
waves (32 < k < 180), corresponding to direct cascade. 
We suppressed condensate by including "artificial" dissi- 
pation on large scales {k < 10). Resulting spectrum is 
given in Fig. [T] (line with short dashes). The best wave 
to see the difference in characteristic waves' scale is to 
have a look at the surface with and without condensate 
(Fig. Surface). Compensated spectrum for direct cascade 
is given in Fig. [3] (right). As one can see exponent of the 
spectrum changed and is now closer to the results of weak 
turbulent theory. The light difference may be a result of 
the influence of the left edge of inverse cascade, which can 
play a role of condensate for short scales corresponding 
to the direct cascade. 

Qualitative explanation of the condensate's influence 
on the short waves could be the following: let us con- 
sider propagating stationary wave with some given slope 
of the front, much longer wave can be treated as a pres- 
ence of a strong background flow. If the direction of the 
flow is opposite to direction of wave propagation the slope 
of the wave's front will increase. This is what we se e in 
our simulations. Average steepness /i = \J < |V?7p > has 
changed: with condensate fic — 0.14, without conden- 
sate fine — 0.12. More detailed picture is given by prob- 
ability distribution functions (PDFs) for surface slopes 
(Fig. [5]l6l) . Although maximums of distributions are well 
described by Gaussian distribution (which is one of the 
assumption of the weak turbulence theory) , we have sig- 
nificant non-Gaussian tails and, what is more important, 
widths of PDFs are different. It means, that in the pres- 
ence of condensate steeper waves are more probable. 
Surface elevation PDFs, given in Fig. [7] in both cases are 
in a good agreement with Tayfun distribution [sot , which 
is the first nonlinear correction to Gaussian distribution. 
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FIG. 3: Compensated direct cascade spectra C(|afe|^)fc° with (left) and without (right) condensate. 
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FIG. 4: Surface of the fluid r\(r) with (left) and without (right) condensate. 



In nature it will result in stronger "whitecapping" : for- 
mation of white foam cap on the crest of the wave caus- 
ing additional transport of energy to the small dissipa- 
tive scale. In the framework of our model such micro- 
wavebreaking is impossible. Dissipation in the system 
prevents formation of strong spectrum tails correspond- 
ing to formation of discontinuities on the surface. Never- 
theless, the mechanism is quite similar: higher steepness 
means stronger nonlinearity in our system. In this case 
for harmonics close enough to the dissipation region gen- 
eration of second and third harmonics acts as fast and 
effective additional process of energy transport to the 
dissipative scales. Processes, corresponding to multiple 
harmonics generation are non-resonant and they are ne- 
glected in the theory of wave turbulence. Also it explains 
why in the experiment in the framework of Zakharov's 
equations ^25] spectra were close to weak turbulent. Za- 
kharov's equations take into account only resonant inter- 
actions and do not describe multiple harmonics genera- 
tion. We can see, that catastrophic events, like formation 
of sharp crests, which cannot be described in the statisti- 
cal framework of kinetic equation, can significantly affect 



physics in the system. The waves kinetic equation can 
be augmented by additional dissipation term to simulate 
this dissipation. As it was shown in recent open field [sij 
and numerical (20| experiments, whitecapping dissipation 
is a phenomenon similar to a second order phase tran- 
sition, so even such a moderate change of the average 
steepness as we observed can cause significant altering of 
the energy transfer mechanism. Our results in the pres- 
ence of condensate can be considered as the first proof 
of a conjecture [3], that Phillips spectrum corresponds 
to a physical picture when a balance between nonlinear 
transport terms and intermittent dissipation takes place. 

-Conclusion - In this Letter author presented results of 
the first direct numerical simulation of the direct cascade 
in the presence of inverse cascade and condensate. The 
importance of condensate as a factor, which increases 
average steepness and stimulates additional intermittant 
dissipation, is demonstrated. Qualitative explanation of 
observed spectra is given. The quantitative explanation 
is a subject of further investigations. Still there is no 
comprehensive theory of whitecapping, which includes 
analysis of fully nonlinear equations. One can use pre- 
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FIG. 5: PDFs of Va;?7 with (left) and without (right) condensate. 
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FIG. 6: PDFs of V^ry with (solid) and without (dashed) con- 
densate. 



sented results for explanation of observed differences in 
the spectra in open see and water tanks experiments. 
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FIG. 7: PDFs of surface elevation rj{f) with (left) and without (right) condensate. 
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